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ABSTRACT 

We present an approximate form for the weak lensing magnification distribution of standard candles, 
valid for all cosmological models, with arbitrary matter distributions, over all redshifts. Our results are 
based on a universal probability distribution function (UPDF), P{rj), for the reduced convergence, rj. 
For a given cosmological model, the magnification probability distribution, P(^), at redshift z is related 
to the UPDF by P(p) = P(rj)/2\K min \, where r\ = 1 + (j« - l)/(2|« min |), and n min (the minimum 
convergence) can be directly computed from the cosmological parameters (fi m and f^)- We show that 
the UPDF can be well approximated by a three-parameter stretched Gaussian distribution, where the 
values of the three parameters depend only on the variance of r\. In short, all possible weak lensing 
probability distributions can be well approximated by a one-parameter family. We establish this family, 
normalizing to the numerical ray-shooting results for a ACDM model by Wambsganss et al. (1997). Each 
alternative cosmological model is then described by a single function £, v (z). We find that this method gives 
P(n) in excellent agreement with numerical ray-tracing and three-dimensional shear matrix calculations, 
and provide numerical fits for three representative models (SCDM, ACDM, and OCDM). Our results 
provide an easy, accurate, and efficient method to calculate the weak lensing magnification distribution 
of standard candles, and should be useful in the analysis of future high-redshift supernova data. 

Subject headings: cosmology: observations — cosmology: theory — gravitational lensing 



1. INTRODUCTION 

The luminosity distance-redshift relations of cosmologi- 
cal standard candles provide a powerful probe of the cos- 
mological parameters Hq, fl m , and (Garnavich et al. 
1998a; Perlmutter et al. 1999; Wang 2000b; Branch et al. 
2001), as well as of the nature of the dark energy (Gar- 
navich et al. 1998b; White 1998; Podariu & Ratra 2000; 
Waga & Frieman 2000; Maor et al. 2001; Podariu, Nu- 
gent, & Ratra 2001; Wang & Garnavich 2001; Wang & 
Lovelace 2001; Kujat et al. 2002). At present, type la su- 
pernovae (SNe la) are our best candidates for cosmological 
standard candles (Phillips 1993; Riess, Press, & Kirshner 
1995). The main systematic uncertainties of SNe la as cos- 
mological standard candles are weak gravitational lensing 
(Kantowski et al. 1995; Frieman 1997; Wambsganss et al. 
1997; Holz 1998; Holz & Wald 1998; Wang 1999; Valageas 
2000a,b; Munshi & Jain 2000; Barber et al. 2000; Premadi 
et al. 2001), and luminosity evolution (Drell, Loredo, & 
Wasserman 2000; Riess et al. 1999; Wang 2000b). Future 
SN surveys (Wang 2000a, SNAP 1 ) could yield thousands 
of SNe la out to redshifts of a few. Since the effect of weak 
lensing increases with redshift, the appropriate modeling 
of the weak lensing of high-redshift SNe la will be impor- 
tant in the correct interpretation of future data. In ad- 
dition, with high statistics it may be possible to directly 
measure the lensing distributions, and thereby infer prop- 
erties of the dark matter (Metcalf & Silk 1999; Seljak & 
Holz 1999). 

In general, determining the magnification distributions 

1 see http://snap.lbl.gov 



of standard candles due to weak lensing is a laborious 
and time-consuming process, involving such techniques as 
ray-tracing through N-body simulations or Monte-Carlo 
approximations to inhomogeneous universes. Here we 
present an easy, accurate, and efficient method to calcu- 
late the weak lensing magnification distribution of stan- 
dard candles, P(n). Our method avails itself of a universal 
probability distribution function (UPDF), P(rj), which we 
fit to a simple analytic form (normalized by the cosmolog- 
ical N-body simulations of Wambsganss et al. (1997)). All 
weak lensing magnification probability distributions, for 
all cosmological models over all redshifts, can then be ap- 
proximated by a one-parameter family of solutions. The 
underlying fundamental parameter is £„, the variance of 
the reduced convergence, rj. To determine the magnifica- 
tion PDF for a given model it is thus sufficient to deter- 
mine £17(2) for that model. We demonstrate this method 
with a number of examples, and provide fitting formulae 
for three fiducial cosmologies (see Table 1). 

Table 1 



Three fiducial models 



Model 


S2m 


n A 


h 


0"8 


SCDM 


1.0 


0.0 


0.5 


0.6 


ACDM 


0.3 


0.7 


0.7 


0.9 


OCDM 


0.3 


0.0 


0.7 


0.85 



2. WEAK LENSING OF POINT SOURCES 

Due to the deflection of light by density fluctuations 
along the line of sight, a source (at redshift z s ) will be 
magnified by a factor fi ~ 1 + 2k (the weak lensing limit), 
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where the convergence k is given by (Bernardeau, Van 
Waerbeke, & Mellier 1997; Kaiser 1998) 
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with 
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and where fi fe = 1 — Q m — fl^, and "sinn" is defined as 
sinh if Qk > 0, and sin if fife < 0. If fife = 0, the sinn and 
fife's disappear. The density contrast S = (p — p)/p. Since 
p > 0, there exists a minimum value of the convergence: 



Or, 



(2) 



The minimum magnification is thus given by p 

1 + 2K m i n . 

Now we define (Valageas 2000a) 

_ Iq 3 dxMx,X, 
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Note that 77 is the average matter density relative to the 
global mean, weighted by the gravitational lensing cross 
section of a unit mass lens along the line of sight to 
the source. This is the same as the direction-dependent 
smoothness parameter introduced by Wang (1999) in the 
weak lensing limit (Wang, in preparation). 

The variance of 77 is given by (Valageas 2000a, b) 
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where A 2 (k,z) — 47rfc 3 P(fc,z), k is the wavenumber, and 
P(k, z) is the matter power spectrum. The window func- 
tion W(Vk9 ) = 2J 1 (Vk6 Q ) / (Vk9 Q ) for smoothing angle 
$o- Here J\ is the Bcsscl function of order 1. Using the hi- 
erarchical ansatz to model non-linear gravitational cluster- 
ing (Balian & Schaeffer 1989), Valageas (2000a, b) showed 
that 



P(V) = f 



27ri£ J) 



,[vy-<j>r,(y)]/tr, 



(5) 



where 4> n (y) — J °° dx (1 — e~ xy ) h{x). The scaling func- 
tion h(x) can be obtained from numerical simulations of 
large scale structure. For iCl, h(x) oc x w ~ 2 (Valageas 
2000a), where lo is the scaling parameter. The uncertainty 
in Eq.(5) comes primarily from the uncertainty in u>. We 
found that the scaling function given by Valageas (2000a) 
leads to large errors in P(i]) for small making it less use- 
ful for calculating P{p) at higher rcdshifts. Although P(p) 
becomes increasingly broad as source redshift increases, 



P(rj) becomes increasingly narrow, since the universe be- 
comes more smoothly distributed at high z (Wang 1999). 
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Fig. 1. — ^J~£^ (for smoothing angle 8 = 1') and — K mi , 
for the three cosmological models of Table 1. 

3. THE UNIVERSAL PROBABILITY DISTRIBUTION 
FUNCTION 

Munshi & Jain (2000) showed that P(rj) is independent 
of the background geometry of the universe, as can be 
seen from equation (5): since weak lensing contributions 
are dominated by a narrow range of the matter power 
spectrum, the scaling function h(x) is independent of cos- 
mological parameters, and hence P(rj) has no explicit de- 
pendence on cosmology (Munshi & Jain 2000). Thus the 
cosmological dependence of P(r)) enters entirely through 
the variance, £ v . We can determine the functional form of 
P(v\€ti) ^ fitting it to accurate calculations of P(p) for 
any cosmological model. The amplification distribution, 
P(p), for arbitrary alternative cosmological models can 
then be found by computing the appropriate K m in [equa- 
tion (2)] and Utilizing p = 1 + 2\K min \(i] — 1) we find 



(6) 



We call P(r]) the universal probability distribution func- 
tion (UPDF), as this one-parameter family of solutions 
underlies all weak lensing magnification PDFs for all cos- 
mologies, at all redshifts. We expect our results to be valid 
in the weak lensing limit, for n < 0.2. In particular, the 
PDFs derived using our formulae are not expected to have 
accurate high magnification tails. 

Figure 1 shows \ f^q and —n m i n computed using equa- 
tions (4) and (2), for the three cosmological models from 
Table 1. For illustration, we give accurate fitting formulae 
in Table 2 for the curves in Figure 1. 



Table 2 

Fitting formulae for curves in Fig.l: 
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SCDM 


.032 


.986 


-.452 


.114 


-.025 


.667 


.482 


-.337 


ACDM 


.021 


1.384 


-.642 


.147 


-.015 


.280 


.766 


-.426 


OCDM 


.032 


1.761 


-.648 


.146 


-.004 


.121 


.703 


-.538 



Universal distribution for weak lensing amplification 



3 



l 

0.8 

^peak 0.6 

0.4 - 




xWambganss et al. 1997 ACDM 
■ •'Munshi & Jain 2000 SCDM, ACDM, OCDM 
• Barber et al. 2000: ACDM 
^Eq.(5) with scaling parameter cj-0.3 ±0.05 
0.2 \~ °OMsothermal Spheres: SCDM, ACDM, OCDM 
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Fig. 2. — The dependence of the three UPDF fitting pa- 
rameters, as functions of y/f^. Data for a variety of different 
models, over a range of redshift, is shown. In addition, best 
fit curves to the Wambsganss et al. points are superposed (see 
equation 8). 



We extract the UPDF, P{ri), from ray-tracing within 
the large scale structure simulations of Wambsganss et al. 
(1997). We then fit the UPDF to the stretched Gaussian 
(Wang 1999): 



p (v\Q = C norm exp 



T] ?]peak 
W T\ q 



(7) 



where C norm , rj pea k, w, and q depend solely on £ v and are 
independent of r\. Note that although eq.(7) accurately de- 
scribes the shape of the UPDF in the range of r] which is 
relevant to weak lensing, one must impose self-consistency 
by restricting 77 < rj max , with rj max chosen such that eq. (7) 
gives the correct Typically, Tjrnax ^ 3-7forl<z<3 
in a ACDM model. We find that (r?) = 1 (hence (pt) = 1) 
for r\ < r\ m ax- This is because rj max is sufficiently large 
so that the contribution of the high 77 tail to the mean is 
negligible. C norm (£, v ) is a normalization constant, chosen 
so that J ?7max P(rj)dri = 1. Figure 2 shows rj pea k, w, and 

q as functions of The points denoted by crosses are 

extracted from the numerical P{(i\z) by Wambsganss et al. 
(1997) for a ACDM model with Vt m = 0.4, n A = 0.6; the 
solid curves arc (x 2 minimizing) best fits to the crosses: 

/ , \ 2 



Vpeakitv) = L002 - L145 - 20 A27 



w(Z v ) = -028 + 3.952 - 1.262 



1 



q(Z v ) = .702 + .509 — — + .008 



, (8) 



The parameter r]p ea k(£,rj) indicates the average smooth- 
ness of a universe; it increases with decreasing £ v (i.e., in- 
creasing z) and approaches r] pea k(^n) = 1 for ^ — > 0. The 



parameter w(^ v ) indicates the width of the distribution in 
the smoothness parameter 77; it decreases with decreasing 
£ v (i.e., increasing z). The £ v dependences of 

(67) and 

w(^ v ) are as expected because as we look back to earlier 
times, lines of sight sample more of the universe, and the 
universe becomes smoother on average. The parameter 
q(£ v ) indicates the deviation of P(r]\t; v ) from Gaussianity 
(which corresponds to q = 0). 

4. COMPARISON WITH OTHER PUBLISHED RESULTS 

The universal probability distribution function, P{rj), 
encapsulated in cqs. (7) and (8), can be used to deter- 
mine the magnification probability distribution, for 
arbitrary cosmological models at arbitrary redshifts. For 
each parameter and redshift, the single free parameter ^ 
determines the full probability distribution. 
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Fig. 3. — The amplification probability distribution, P(/J,), 
derived from ray-tracing simulations by Munshi & Jain (2000) 
(open symbols) for smoothing angle 60 = 1', source redshift 
z s = 1, and the three cosmological models of Table 1, together 
with P{n) computed using our UPDF, with K miTl and £ v com- 
puted using equations (2) and (4). 
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Fig. 4. — The P(n) from three-dimensional shear matrix 
calculations of N-body simulations by Barber et al. (2000) 
(circles) for a ACDM model with 57 m = 0.3, £Ia = 0.7 at source 
redshifts z s — 3.6, 2, 1 (peaking from left to right), together 
with P(fi) computed using our UPDF, with Kmin computed 
using equation (2) and l; v inferred from Table 4 of Barber et 
al. (2000). 

Figure 3 shows the P{n) from ray- tracing simulations by 
Munshi & Jain (2000) for smoothing angle 6 = 1', source 
redshift z s = 1, and three cosmological models from Table 
1, together with P(/i) computed using our UPDF for the 
K m i n and ^ v computed using equations (2) and (4). 

Figure 4 shows the P(fi) from three-dimensional shear 
matrix calculations of N-body simulations by Barber et 
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al. (2000) for a ACDM model with tt m = 0.3, Q A = 0.7 
at source redshifts z s = 1, 2, 3.6, together with P(n) 
computed using our UPDF, with K m i n computed using 
equation (2) and ^ inferred from Table 4 of Barber et 
al. (2000) [equation (4) was not used to compute £ v due 
to our lack of knowledge of the smoothing angle 6 that 
corresponds to their results]. 

Our UPDF gives P{p) in excellent agreement with N- 
body calculations. To make this more apparent, we have 
extracted P(rj) from the P(/x) obtained via N-body calcu- 
lations by Munshi & Jain (2000) and Barber et al. (2000), 
and fitted them to the functional form of equation (7) . The 
resultant coefficients are plotted in Figure 2. There is very 
good agreement in the peak location ?7 pe afc(C»?) an d width 
indicator w^), but larger scatter in the non-Gaussianity 
indicator q(£ v ) extracted from Munshi & Jain (2000) and 
Barber et al. (2000). The latter may arise partly due 
to the fact that in both cases we poorly resolve the non- 
Gaussian tails, which are crucial to determining accurate 
values of q. In addition, the weak lensing condition breaks 
down for the high /x tails, which could be significant for 
small £ v . Also plotted in Figure 2 are the coefficients ex- 
tracted from fitting the analytically computed P(rj) [see 
Eq.(5)], following Munshi & Jain (2000), for the scaling 
parameter lo — 0.3 ± 0.05. These P{rj) have not been 
tested for z > 1 (i.e., for small although the devia- 
tions are expected to be small, since P(rj) peaks close to 
r\ = 1 at small £ v [see equation (7)]. 

Figure 2 also shows the P{rj) coefficients extracted from 
ray-tracing of randomly placed singular isothermal spheres 
(SIS), following the prescription of Holz & Wald (1998), for 
the three cosmological models of Table 1; these are in good 
agreement with the fitted coefficients from Wambsganss et 
al. (1997). We find that improved statistics leads to bet- 
ter agreement between the q(^ v ) from our ray-tracing of 
randomly placed mass distributions and that from Wamb- 
sganss et al. (1997), while having much less impact on 
Vpeak{^7j) and w(£ v ). This is as expected, since improved 
statistics fills out the non-Gaussian tails of the P(fjb). 

5. SUMMARY AND DISCUSSION 

We have derived a simple and accurate method to com- 
pute the weak lensing magnification distribution, P(n), for 



standard candles placed at any redshift in arbitrary cosmo- 
logical models. We use a universal probability distribution 
function (UPDF), P{rj\£, v ), which is independent of cosmo- 
logical model; the dependence on cosmology entering only 
through the variance, of the reduced convergence, r\. 
The UPDF is fit accurately by a 3-parameter stretched 
Gaussian distribution [cq. (7)]. We give polynomial fitting 
formulae [eq. (8)] for the three parameters rj pea k(^ v ) (aver- 
age smoothness), w(^) (smoothness variation), and q(^n) 
(non-Gaussianity), which we normalize to the N-body sim- 
ulations of Wambsganss et al. (1997). The magnification 
PDF, P(fi,z), can then be determined from the UPDF 
using equation (6). We expect our results to be valid in 
the weak lensing limit, for k < 0.2. The extension of our 
method to high magnifications will be presented elsewhere. 

To test the robustness of this method, we have compared 
our results against three alternate independent methods 
(see Fig. 2). We find excellent agreement with: (1) the 
N-body calculations of Munshi & Jain (2000) and Barber 
et al. (2000), with some scatter in the non-Gaussianity 
indicator </(£,,), which is consistent with the limited statis- 
tics at low z (i.e., large £,,) of these N-body results and the 
breaking down of the weak lensing condition at high /x, (2) 
the analytical calculation [see equation (5)] following Mun- 
shi & Jain (2000), where the latter has been verified by ray 
tracing experiments, and (3) the ray-shooting of randomly 
placed SIS mass distributions, following the prescription of 
Holz & Wald (1998). 

We expect these simple, universal forms for the weak 
lensing distribution to be useful in addressing high red- 
shift data, and in particular, in the analysis of results from 
future supernova surveys. 
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